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Abstract 

We consider a preferential growth model where particles are added one by one to 
the system consisting of clusters of particles. A new particle can either form a new 
cluster (with probability q) or join an already existing cluster with a probability 
proportional to the size thereof. We calculate exactly the probability Vi(k,t) that 
the size of the i-th cluster at time t is k. We applied our model as a background for 
a microscopic economic model. 



1 Introduction 



Preferential growth describes systems consisting of groups of entities where 
the probability of the attachment of a new entity to one of the groups is an 
increasing function of the group's size. Recently these types of models were 
used to describe networks like the World Wide Web (WWW) [1], Internet [2], 
statistics of scientific citation [3] etc which seems to have the common property 
of scale invariance. It turned out that this behavior is due to two main facts 
namely that they are continuously evolving and that the attachment probabil- 
ity is proportional to the group's size. Many features of the above mentioned 
systems were analyzed and also analytic calculations have been presented for 
the most important quantities [4] but they are only valid in the asymptotic 
time case. Here we present the full-time dependent solutions. 

We applied this kind of model to describe herding in an economic system. It 
was suggested [5] that large price movements on the stock market can be ex- 
plained by the fact that traders are not individual participants but they prefer 
to form groups, within a group each trader share the same strategy. Our as- 
sumption is that those groups are created according to a preferential growth 
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model. Preferential growth occurs naturally in economics [6]. The power-law 
behavior of the resulting group distribution can account for the fat-tail dis- 
tribution of the return which is one of the stylized fact that characterize the 
stock market. 

In the paper we present the full-time dependent solution of the growth model, 
and briefly describe the most important results of our economic model. 



2 Growth model 

The system consists of individual groups of entities. The initial condition is 
one group with one entity in it. At each time step one new entity is added 
to the system. With probability q it creates a new group, with p — 1 — q it 
will belong to one of the existing groups. The probability that it joins the ith 
group is proportional to the ith group's size (ki/N). The system's dynamic 
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Fig. 1. Illustration of the dynamic of our model. 

- the time evolution of the ith group's size distribution, Vi(k,t) - can be 
described by a master equation: 

Vi(k,t) = p^^-V t (k - l,t - 1) +p (l - v ^ t ~ ^ + 

+ (1 - p) Vi(k, t-l) + (l-p) lU-^t - 1) " hi)- (!) 

The quantity Hi(t) denotes the probability that at time t the system consists 
of i group: n i (t) = (^Z]) p 1 ^ 1 ^^ 1 ^! —p) l ~ l . Our goal is to determine the 
analytic full-time dependent form of the following quantities: size distribution, 
Vi(k, t), and mean, (ki)(t), of the individual groups and the average group size 

1 * 

distribution, P(k, t) = - J~]Vi(k, t). Here we do not go into the details of the 

* i=i 

deduction we only present the results. (For more informations see [7]). For the 
size distribution of the individual groups one gets: 
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From Vi(k,t) we can easily calculate the analytic form of P(k,t) and (ki)(t) 
by simple substitution in their definitions and we get: 
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In the case when the system is large (t — > 00) it is worth to analyze the 
asymptotic case of the above quantities in order to gain a much simpler form. 

i) Individual group size distribution 

F(t - lp) 

In Eq. (2) there are two time dependent terms. One is — - — —, the 

r(t)r(l - lp) 

other is the sum in the bracket. For large t values the second term will be 
proportional to a hypergeometric function, 2F±(i,i — l;i — lp;p), and will be 
time independent. The first term will converge to t~ lp which is a fast decaying 
function of I so for t^>k one can assume that only the first term of the sum 
in Eq. (2) gives non-negligible component: 



hm7MM)=t- p (l-P) i_1 p^T 

t-KX> Y(l — p) 



2 F l (t,t-l;t-p;p) + 0(t- 2p ) (5) 



For large i values the above form simplifies further and we finally get: 



lim Vi(h,t) = , , 



(6) 



ii) Mean of individual group size 



In the analysis of the mean value we assume that for k <^ t the distribution 
can be described by Eq. (5) for larger k values it has a fast decay so: (h)(t) 
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k Vi(l,t — > oo). (The value k* can be denned e.g. as the inflection point 

k=i 

of Vi(k,t)). Here again analyzing the large i case we get: 
iiij Distribution of average group size 

In the determination of the asymptotic case of the average group size distri- 
bution, P(k), we can apply two ways. One is to simply analyze the limit of 
Eq. (3) or we can also start directly from the master equation, Eq. (1), by 
summing it up for i — 1 ... t. In both cases we get that in the long time limit 
the distribution will be time independent and will have a power-law decay: 

V{k ) = -j^ £fi— ? ~ k /p . (7) 

Here we would like to mention that however our growth model is not a network 
model for particular values of the parameter (p = 0.5) it can be considered as 
a mean field analogy of the Barabasi's network model. The resulting exponent 
of the distribution, (—3), agrees with their result [1]. 



3 Economic model 

The main properties of stock markets price movements are characterized by 
the so called stylized facts, the distribution of the logarithmic return fat-tailed, 
autocorrelation of the return is short range while that of the absolute value 
of return decays with a power-law. We try to reproduce these properties by a 
simple model. 

Power-law decay of the distribution of the logarithmic return can be explained 
by the herding behavior on the market [5]. Price changes are assumed to be 
proportional to the excess demand - the difference between the amount of 
buying and selling orders. If traders form groups so that within each group ev- 
eryone shares the same strategy than the logarithmic return can be expressed 

n 

as: x ~ ^Sj0j, where 0« is the strategy (to buy, to sell or not to trade), Sj 
i=i 

is the size of the ith group. If the number of groups that are active (trading) 
at a given time step is small the distribution of the logarithmic return, P(x) 
will be similar to that of the groups' size distribution, P(s). In our model 
groups are created according to the above mentioned growth model. Due to 
Eq. (7) P(s) will have a power-law decay which indicates that for law activity 
the distribution of logarithmic return is also a power-law. Similar results were 
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get by [8] but they considered trading groups as two dimensional percolation 
clusters. 

Power-law correlation of the absolute value of price is due to the long range cor- 
relation of the trading volume (number of trades in a given time interval) [9], 
which is in our model the sum of the active groups' sizes, S a = Yn=i s i- We 
defined the activity (fraction of groups that are trading at a given time step, 
n a /n) as the function of the deviation of the actual price from some funda- 
mental price, p ; see Fig. 2. (We took the fundamental price for constant over 
time.) We ensure with this definition that the activities of the proceeding 
time steps are correlated. For a given activity value the corresponding trading 
volumes can be different but in average larger activity implies larger trading 
volume which indicates that the sum of the trading groups' sizes will also be 
correlated. The form of the function was chosen intuitively, it is based on the 
hypothesis that small deviation from the fundamental price inspire traders - 
their activity rises - while for large deviations traders become careful so their 
activity falls back. 
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Fig. 2. Price dependence of the activities. a + is the fraction of the buying groups, 
a_ is the fraction of selling groups at a given time step. For ln[p(t)/po] < the 
activities a + and a_ are equal. If ln[p(t)/po] > b* then a+ = 0, a_ = ao while for 
ln\p(t)/po] < b* => a_ = 0, a + = ao- ao is the reciprocal of the number of groups. 

The results of our model can be seen in Fig. 3. It is surprising that for the 
distribution of the logarithmic return we not only got a simple power-law but 
a distribution with two different exponents (Fig. 3/a) which is similar to what 
was measured on the real market [10]. This is, however, only a qualitative 
agreement, the value of the exponents differ from the measured ones [10]. 
To analyze the autocorrelation of the absolute value of return we study the 
variance of the sum of these variables (normalized to unit variance): 



n-l 



D 2 (Y, x i) = D 2 (n) = n + 2 Y.in - m) C{m) ~ n 2 \ 



i=l m=l 

where C(m) ~ mr K is the autocorrelation function. For k > 1 we have 5 = 0.5, 
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Fig. 3. a) Distribution of the logarithmic return. It has on the log-log plot two 
approximately linear regimes. The first has an exponent ~ 1.72 the second ~ 3.36 b) 
Square root of the variance, D(n). The dashed line corresponds to the uncorrelated, 
5 = 0.5, time series. The full line is a fit with exponent 5 = 0.95. 

for variables with long-range correlation (k < 1) 0.5 < 5 < 1, and the corre- 
lation exponent, k, can be calculated through k = 2 — 25; see Fig. 3/b. It is 
clearly shown that the time series of the absolute value of the logarithmic re- 
turn is correlated but the exponent is much smaller (~ 0.05) than the reported 
one (~ 0.3) [10]. 

The advantage of our model is its simplicity and that the results have a qualita- 
tive agreement with the stylized facts. However, there are tunable parameters 
in it and the results - the values of the exponents of the distribution and of 
the autocorrelation - depend on them. Further study is needed to clarify the 
role of these parameters. 
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